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1. Introduction 

We study here the convergence of solutions of the partial differential equation 
(e _» o+, 5 = 5(e) -» 0) 

u t + f(u) x = eu xx + 5u xxx , u = u € (x,t), x £ R, t > 0, (1.1) 

toward weak solutions of the corresponding hyperbolic conservation laws: 

u t + /(u)a; = 0, u = u(x, t), X £ R, t > 0, (1.2) 

where the flux / : R — > i? is a smooth function with (at most) linear growth at 
infinity, that is, for some M > 

\f'(u)\<M, u£R. 

Equations of the form (1.1)-(1.2) arise in fluid dynamics when both viscosity 
(diffusion) and capillarity (dispersion) play a role. The diffusion e smoothes out 
the discontinuous solutions of (1.2), while the dispersion 5 causes high-frequency 
oscillations. 

In this paper, we establish that the solutions u e of (1.1) converge toward a weak 
solution of (1.2) provided 

5 = 0(6 2 ). (1.3) 

When the stronger condition 

5 = o(e 2 ) (1.4) 

holds, we prove that the limit coincides with the entropy solution determined by 
Kruzkov's theory [10]. We point out that these conditions are sharp since, in the 
limiting case, 

5 = Ke 2 for some K G R, (1.5) 

limiting solutions may violate Kruzkov's entropy conditions [8, 5, 13, 2]. Further- 
more, when (1.3) is violated, the solutions are highly-oscillatory and fail to converge 
in any strong topology as noted by Lax and Levermore [12]. (See also Lax [11]. ) 

The singular limit problem above was first tackled by Schonbek [19], who estab- 
lished the optimal rate (1.3) for Burgers equation, that is, 



u 2 



/(«) 

and for the class of flux-functions 

u 2p+l 

f(u) = - , P >1. 

Jy 1 2p + l 

She also gave a convergence result for general fluxes with quadratic growth at in- 
finity, however under the stronger condition on 8 = 0(e 3 ). As another important 
contribution in [19] , Schonbek introduces a generalization of the method of compen- 
sated compactness (Tartar [21] and Murat [18]) allowing to handle sequences that 
are bounded in LP for finite p > 1 only. Next, following [19], LeFloch and Natalini 
[15] studied equations like (1.1) but with nonlinear (even singular) diffusion, and 
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established strong convergence results toward entropy solutions of (1.2). See also a 
convergence result for systems in Hayes and LeFloch [6]. 

In the second part of this paper, we also deal with the convergence of solutions 
of multi-dimensional equations similar to (1.1)-(1.2). For multi-dimensional equa- 
tions, the compensated compactness method no longer applies and the proofs are 
based instead on DiPerna's uniqueness theory for entropy measure- valued solutions 
(DiPerna [4], Szepessy [20], and Kondo and LeFloch [9]). Our approach is similar to 
Correia and LeFloch [3] where nonlinear diffusion terms are treated under a strong 
assumption on the ratio of the dispersion to the diffusion. 

To summarize, the main contribution in the present paper is the derivation of 
a priori estimates (Theorems 2.1 and 3.1) which cover general flux-functions (with 
at most linear growth at infinity) and lead to an optimal condition on the balance 
between the diffusion and the dispersion (Theorems 2.2 and 3.2). 

Further material on classical and nonclassical entropy solutions generated by 
diffusive-dispersive limits can be found in [1, 2, 5, 6, 7, 8, 13, 14, 15, 16, 17, 19] 

2. One-Dimensional Conservation Laws 
Consider a family u e of smooth solutions to 

u t + f{u) x = eu xx + Su xxx , u = u £ (x,t), (2.1) 

u(x,0) =u e (x), xeR, (2.2) 

where e — > 0+ and S = 5(e) — > 0. Under suitable conditions on the initial data 
Uq : R — > R, the solutions (and their derivatives) decay rapidly at infinity, so that 
all the a priori estimates given below are rigorously justified. We want to show that 
the solution of (2.1)-(2.2) converges toward a weak solution of the problem 

«t + /(«)* = 0, u = u e (x,t), (2.3) 

u(x,0) = u (x), x£R, (2.4) 

where no : R — > R is a given initial data. A minimum requirement is the weak 
convergence (for instance in L 2 (R)) 

which is always assumed throughout this paper. The following convergence theorem 
covers both cases where the diffusion are in balance or dominates the dispersion. 

Theorem 2.1. Suppose that the flux-function the flux-function f is Lipschitz con- 
tinuous on R and that the initial data u belong to L 2 (R) . Then the solution u e of 
(2.1)-(2.2) satisfies the following a priori estimates: 

\\u e (t)\\L2{R) < \\u € \\ L 2{R), t > 0, (2.5a) 
oo ||«olk 2 (H) + v 7 ^ IKJIl 2 cr), t > 0, (2.5c) 
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and 

^\\ U ^\\ L 1^ R + >L 2 (R ^ ^ V 2 ll/'lloo ||t*0lU a (H) + ^ ll«LIU 2 (J?)- ( 2 - 5d ) 

Proof. Throughout the calculation and for simplicity, we omit the upper-index e. 
To any smooth function U : R — > R we can associate a "flux" i* 1 : i? — > R by 
F» = U'(u)f'(u), u£R. Multiplying (2.1) by £/» we find 

U(u) t + = e (tf'faK)* - e tf"(u) + 5 m^)^ - «* W 

Integrating over the whole space, it follows that 

^ / £/(u)cte + e / U"(u)u 2 x dx = 5 [ U"(u)(^] dx 

dt JR JR JR V 2 / x ^ 6) 

= - S -J R U"'(u)uldx. 
Integrating in time over some interval (0,i), we arrive at the general identity: 
U(u(t))dx + e [ [ U"{u)u 2 x dxdt= [ U(u )dx-^ [ [ U"\u)u 3 x dxdt. 

JO JR JR 2 JO JR 



R 



(2.7) 

Choosing first U{u) = u 2 in (2.7), we see that 

f u(t) 2 dx + 2e [ [ u 2 x dx= f u 2 (x), (2.8) 

JR JO JR JR 

which gives immediately (2.5a) and (2.5b). 

Next, we differentiate (2.1) with respect to x and we multiply by u x : 

- (u 2 x ) t + u x (f'(u) u x ) x =e (u x u xx ) x -eu 2 xx + 5 (u x u xxx - -u 2 xx ) x . 

Integrating in space, we get 
1 d 



2 dt 



I uldx + e I u 2 xx dx= [ u xx f (u) u x dx = ~ [ f"(u)u x dx. 

JR. JR JR * JR 

Hence, integrating over some interval (0, t), we find 

f u x (t) 2 dx + 2e [ [ u 2 xx dxdt= [ u 2 0x dx- [ [ f" {u)u 3 x dxdt. (2.9) 

JR JO JR JR JO JR 

Multiply (2.9) by 5 and add it up with (2.7): 
5 / u x (t) 2 dx + 2eS / / u 2 x dxdt 

JR JO JR 

= [ U(u )dx- [ U(u(t))dx + 5 f u 2 Qx dx-e f [ U"{u)u 2 x dxdt 

JR JR JR JO JR 

-5 f [ f"{u)uldxdt- 5 - [ [ U"'{u)u 3 x dxdt. 

JO JR 2 JO JR 
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Choosing U given by 

U(u) = -2 [ U (f(v)-f(0))dv (2.10) 
J o 

the last two terms in the above identity cancel out. Since 

U"(u) 

-c< —±- L < c:= H/'lloo, foralluGi?, 

thus 

-c« 2 < U(u) < cu 2 , for all u & R, 

we finally obtain 

<5 / u x (t) 2 dx + 2 e 5 / / u 2 xx dxdt 
Jr Jo Jr 

< / cu^dx + / cu(t) 2 dx + 5 / UQ a ,(ix + 2e / / cu 2 dxdt. 
Jr Jr Jr Jo Jr 

Hence using (2.8) 

5 / u x (t) 2 dx + 2e5 / / u 2 x dxdt <2c / Uq + 5 / ul x dx, 
Jr Jo Jr Jr Jr 

which leads to (2.5c) and (2.5d). The proof of Theorem 2.1 is completed. □ 

Recall that by Kruzkov' theory, given uq € L 2 (R) the Cauchy problem (2.3)-(2.4) 
admits a unique entropy solution u € L°° (R + ,L 2 {R)) in the sense of Kruzkov's 
theory. See [10, 4, 20, 9]. 

Theorem 2.2. Assume that, for some constant C > independent of e, 

IKIU^K) + v 7 ^ < Co- (2.11) 



(1) ^4s e — > 5 = 0(e 2 ) (a subsequence of) the solution u e of (2.1)-(2.2) 
converges in Lf oc L^ oc (i?)) f/or a// 1 < p < oo and 1 < q < 2) toward 
a weak solution of the problem (2.3) -(2.4). 

(2) If the stronger condition 5 = o(e 2 ) holds, then the limit is the unique entropy 
solution in the sense of Kruzkov. 

In Case (1) a subsequence of u e (at least) converges strongly, while in Case 
(2) the whole sequence converges strongly. We can conjecture that, in fact, the 
whole sequence should converge in Case (1) as well, but proving it would be very 
challenging since it requires a uniqueness result of nonclassical entropy solutions. 
(See also LeFloch [14].) 

Proof. We will apply the general convergence framework established by Schonbek 
[19]. Based on (2.11) and the uniform estimate (2.5a) derived earlier, we can select 
a subsequence of u e converging "in the sense" of the Young measures. To apply 
[19], we only need to control the entropy dissipation measures associated with the 
equation (2.1). Let U be a smooth function with (at most) linear growth at infinity 
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and, more precisely, such that U' and U" are uniformly bounded on R. Consider 
the distribution 

T e = U(u £ ) t + F(u%, 
where as usual F' = U' f. With obvious notation consider the decomposition 

r« = e (IT (« e ) u%) x - e tf"(u e ) (u%f + S (U'(u^ u% x ) x - 5 U"{u*) u% u% x 
= T\+V\+T\+T%. 

The estimates below hold for all smooth function 9 : R x R + — > with compact 
support in (x, i). 

Consider first the term rf . By Cauchy-Schwarz inequality, we get 



/7 r > 

Jo Jr 



9 dxdt 



'o jr 
< eC||i£ 



f eU'(u e ) 
Jr 



u% 9 X dxdt 



x ''l 1 (r + ,l 2 (r)) 11 x "l°°(r + ,l 2 {r)) 

< a -> o, 



(2.12i) 



where we used (2.5b). This proves that Tf converges to zero in the sense of distri- 
butions. 

Next we simply point out that, by (2.5b) again, the second term Y\ remains 
uniformly bounded in L 1 : 



/o jr 

To estimate T 3 we use (2.5d): 



\T\\dxdt < - \\u \\ 2 lHr) . 



(2.12ii) 



n 

Jo Jr 



T% 9 dxdt 



o Jr. 



< &c\Wlx\ 



U'(u £ )u e xx 8 X dxdt 
li(r+,l->(r)) W ex h°°(R + ,L*(R)) (2.12iii) 



0. 



provided that the mild condition 5 = o(e) holds. Therefore T| tends to zero in the 
sense of distributions. 

Finally, we deal with the last term as follows: 



n 

Jo Jr 



H 9 dxdt 



5U"{u e )u% u% x 9dxdt 



'o Jr 

<SC\\uiJ 



01 



l°°(r + ,l 2 {r)) 11 x"L°° (r+,l 2 (r)) 11 11 L°°(RxR + ) 



(2.12iv) 

where we use (2.5b) and (2.5d). The upper bound above tends to zero iff S = o(e 2 ), 
in which case we can conclude that T| tends to zero in the sense of distributions. 
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Under the weaker assumption 5 = 0(e 2 ), we see that T\ is solely bounded in 
L 1 (RxR + ) as is T e 2 . 

The conclusion (1) of the theorem follows immediately from the uniform bounds 
(2.12) by applying Schonbek's convergence theory. Her arguments only show that 
a subsequence of u e converges and that the limit is a weak solution of (2.3)-(2.4). 
On the other hand, assuming now the stronger condition 5 = o(e 2 ) and restricting 
attention to convex functions U, in view of (2.12) again and the expression of T| 
we see that the entropy dissipation decomposes in the form 

where f e — ► in the sense of distributions and T e is a non-positive bounded measure. 
This shows that all of the entropy inequalities hold in the limit e — > 0. Thus the 
limit coincides with the unique entropy solution of the problem. □ 

3. Multi-Dimensional Conservation Laws 

The estimates and the technique of proof in Section 2 do not apply to multi- 
dimensional equations, and markedly different arguments are discussed now. Con- 
sider the following Cauchy problem: 

d d d 

u t +^2fj(u) Xj = e^2u XjXj +S^2u XjXjXj , u = u e (x,t), x G R d , t > 0, (3.1) 

j=i j=i j=i 

u(x,0) =u e (x), x G R d . (3.2) 

Provided the initial data Uq converge weakly to some limit Uq (in L 2 , say), we will 
now prove that the solutions of (3.1)-(3.2) converge toward the entropy solution of 
the associated hyperbolic problem: 

d 

u t + J2fj( u U =°i u=(x,t),xeR d ,t>0, (3.3) 
i=i 

u(x,0) =u Q (x), xeR d . (3.4) 

Precisely our result are as follows: 

Theorem 3.1. Suppose that the flux-function f is Lipschitz continuous on R and 
that the initial datauo belong to L 2 {R d ). Then the solution u e o/(3.1)-(3.2) satisfies 
the following a priori estimates: 

IKWIIl^) < IKHl^); *>0> (3.5a) 

^IKI'l^+.l^)) ^ KIIl^), ( 3 - 5b ) 

for all j = 1, . . . , d and t > 

e|K 3 .(*)llL2(jjd) < v^H/jlloo IKIIi^iid) +e\\u £ o Xj \\mRi) (3-5c) 
and for all j,k = 1, . . . , d 

e3/2 IK J a;JI L1 ( J? + iL 2 (K d ) ) < V^H/jHoo |Klli 2 (K d ) + 6 \\ U lx 3 \\^(R d ) ■ (3.5d) 

For each uq G L 2 (R d ), the Cauchy problem (3.3)-(3.4) admits a unique entropy 
solution « G L°°(R + ,L 2 (R d ) in the sense of Kruzkov. See again [10, 4, 20, 9, 14]. 
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Theorem 3.2. Assume that, for some constant Co > independent of e, 

d 

IKIIl^) Yl IKsJl^h) ^ Co- (3-6) 

Then, when e — > 0+ wit/i 5 = o(e 2 ) ; i/ie solution u e of (3.1)-(3.2) converges in 
Lf oc (R + , Ll oc (R )) (/or all 1 < p < oo and 1 < q < 2) toward the unique entropy 
solution in the sense of Kruzkov of the Cauchy problem (3.3) -(3.4). 

Recall again that the condition 5 = o(e 2 ) is sharp since, in the opposite case, 
nonclassical solutions violating the Kruzkov entropy inequalities could arise in the 
limit. 

Proof of Theorem 3.1. We omit the upper- index e in the following calculation. To 
derive the L 2 bound (3.5a), we multiply the equation (3.1) by u and get 

(^) t +E*i(«)., 

d d ^. d d 

where Fj = u /j is normalized by the condition i*j(0) = 0, j = 1, . . . , d. Integrating 
over space, we get 



d [ , l2 , 

— / |«| dx = —2e 

at J R d 



,. d 



and so for all i > 

/ |u(t)| 2 (ix + 2e / / V |tt Xj . | 2 dxdt = f \u \ 2 dx. (3.7) 

^i? d ^0 JR d j =1 3 JR d 

To estimate the gradient of u, for k = 1, ■ ■ ■ ,d we differentiate the equation 
(3.1) with respect to the variable x^ and then multiply by u Xk . The right-hand side 
of (3.1) is linear in u thus the calculation for this side is identical to the one we 
just made, but with u replaced with u Xk . On the other hand, the flux term in the 
left-hand side is nonlinear and requires a specific argument: 

— \u Xk \ 2 dx-J2 2u XkXj f!>(u)u Xk dx = -2e V / \u XjXk \ 2 dx, 
at J R d j Rd j R d 

so after integration in time 

n ft p 

j \u Xk (t)\ 2 dx + 2e V] / / \u XjXk \ 2 dxdt 

JR d - =1 ■/ 

|n 0a; J 2 (ix + 2||/fc|| oo V] / / |« Xj . Xfc | |u x J dxdt 
• =1 ■/ o .AR d 

/• II f II 2 ^ /"* T d r-t [■ 

< \u 0xk \ 2 dx+ l -^---dy j / KJ 2 dxdt + e V / / |^. x J 2 dx^. 

e =1 >/i? d =1 JO JR d 
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Observe that the last term of the right-hand side coincides with the last term of 
the left-hand side. Therefore, multiplying the above inequality by e 2 and using the 
entropy dissipation bound in (3.7), we deduce that 



/ e 2 \u Xk (t)\ 2 dx + V] / / e 3 \u XjXk \ 2 dxdt 
jR d - =1 J o JR d 

< / e 2 \u 0xk \ 2 dx+ 11/fcH 2 ^ / / de\u Xk \ 2 dxdt 

JR d JO JR d 

< / e 2 Wo Xk \ 2 dx + dWflWlo \ \u \ 2 dx. 

JR d JR d 



(3.8) 



□ 



Proof of Theorem 3.2. We will rely on the convergence framework proposed by 
DiPerna [4] for L°° solutions and generalized to LP solutions by Szepessy [20] and 
Kondo and LeFloch in [9]. 

Consider a Young measure v associated with the sequence u e and based on the 
uniform L 2 bound (3.5a). (Such Young measures are described in Schonbek [19]). 
To show that v is an entropy measure-valued solution, we must check entropy 
inequalities associated with the equation (3.3), that is, 

d 

J=l 

where U : R — ► R is a convex function with (at most )linear growth at infinity and 
the entropy flux Fj = U' /j is normalized so that Fj(0) = 0. 

By the definition of the Young measure, we only need to establish that, in the 
decomposition 

d 

i=i 

d 

= £ d 3 (e U'(u<) djif + 5(e) U'(u<) d 2 ^) 

3 = 1 ^ ' 



we have 
and 



d 

- e U"{u e ) \d jU e \ 2 + 5(e) U"(u e ) d jU e d 

3 = 1 

-n + ri + ri + n, 

rs,ri,n->o 
r^<o. 



These convergence properties precisely were established in the proof of Theorem 
2.2, at least for one-dimensional equations. The extension to multi-dimensional 
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equations is immediate in view of the uniform estimates (3.5). A detained dis- 
cussion of the initial condition at t = (which is based on using suitable entropy 
inequalities) can be found in Kondo and LeFloch in [9]. This completes the proof 
that the convergence framework in [9] applies and provides the strong convergence 
toward the unique entropy solution of (3.3)-(3.4). □ 
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